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ABSTRACT 

We present numerical simulations of driven magnetohydrodynamic (MHD) turbulence with 
weak/moderate imposed magnetic fields. The main goal is to clarify dynamics of magnetic field growth. 
We also investigate the effects of the imposed magnetic fields on the MHD turbulence, including, as a 
limit, the case of zero external field. Our findings are as follows. First, when we start off simulations with 
weak mean magnetic field only (or with small scale random field with zero imposed field), we observe 
that there is a stage at which magnetic energy density grows linearly with time. Runs with different 
numerical resolutions and/or different simulation parameters show consistent results for the growth rate 
at the linear stage. Second, we find that, when the strength of the external field increases, the equilib- 
rium kinetic energy density drops by roughly the product of the rms velocity and the strength of the 
external field. The equilibrium magnetic energy density rises by roughly the same amount. Third, when 
the external magnetic field is not very strong (say, less than ~ 0.2 times the rms velocity when measured 
in the units of Alfven speed), the turbulence at large scales remains statistically isotropic, i.e. there is no 
apparent global anisotropy of order B /v. We discuss implications of our results on astrophysical fluids. 

Subject headings: ISM:general - intergalactic medium - MHD - turbulence 



1. INTRODUCTION 

Most astrophysical fluids arc magnetized. Magnetic field 
in an astrophysical system can be divided into two compo- 
nents: large-scale regular and small-scale random compo- 
nents. The generation of magnetic field may involve with 
two separate issues: generation of the large-scale regular 
field (or mean field) and generation of small-scale random 
field. 

The generation or growth of large-scale regular fields is 
an important topic in astrophysics. However, in this pa- 
per we assume fixed large-scale regular fields are already 
present and we only investigate how small-scale random 
fields are generated from the imposed large-scale fields. 
That is, we investigate the growth of magnetic energy in 
the presence of fixed mean fields. Therefore, the genera- 
tion/growth of large-scale regular field itself is not a topic 
of this paper. Readers interested in the topic may refer 
to mean-field dynamo theories (see Moffatt 1978; Parker 
1979; Krause & Radler 1980; Brandenburg & Subramanian 
2005). The type of magnetic energy growth we deal with in 
this paper is so-called small-scale turbulence dynam^, the 
origin of which can be traced back to the papers by Batche- 
lor (1950) and Kazantsev (1968). We also study how prop- 
erties of magnetohydrodynamic (MHD) turbulence (e.g. 
energy densities, power spectra, global anisotropy) change 
as the strength of the mean field changes. 

When we introduce a mean field in a turbulent medium, 
the mean field interacts with turbulent motions. There 
are two distinct types of interaction, depending on the 



strength of the mean magnetic field. When the large-scale 
regular magnetic field is weak, turbulent motions stretch 
the magnetic field lines and, as a result, the magnetic en- 
ergy density increases. MHD turbulence near the scale of 
the largest energy-containing eddies (i.e, the outer scale or 
the energy injection scale) will be more or less like ordi- 
nary hydrodynamic turbulence with small magnetic back 
reaction. 

MHD turbulence in intracluster medium and intergalac- 
tic medium may fall in this type of turbulence. The origin 
of the seed magnetic fields is still uncertain. However, 
whatever the origin is, turbulence motions can produce 
small-scale field through the stretching of the seed fields. 
In this regime, the growth timescale will be of great im- 
portance. If the growth time scale is shorter than the age 
of the universe in a system, we expect that the system is 
strongly magnetized in the present time. Otherwise, we 
expect that the system is weakly magnetized. 

On the other hand, when the imposed mean magnetic 
field is strong in the sense that the turbulent eddy turnover 
rate at the large scale (i.e. L/v) is slower than the Alfvenic 
rate of the same scale (i.e. £/-BojE the resulting turbu- 
lence can be described through the nonlinear interaction of 
waves. The turbulence can be either strong, meaning that 
the cascading happens within one eddy turnover time or 
weak, meaning that the cascading takes more than one 
eddy turnover time. A classical study of Iroshnikov (1963) 
and Kraichnan (1965) presents an example of weak MHD 
turbulence. This is a hypothetical isotropic MHD turbu- 
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lence, while we know by now that the actual MHD tur- 
bulence is anisotropic (Shebalin et al 1983; Higdon 1986; 
Goldreich & Sridhar 1995; Matthaeus et al. 1998; see also 
Biskamp 2003 and references therein). An exact analyt- 
ical treatment of weak MHD turbulence can be found in 
Galtier et al. (2000). A successful model of the strong 
MHD turbulence was presented in Goldreich & Sridhar 
(1995). 

In this paper we deal with the former case: turbu- 
lence with weak imposed magnetic field. In the interstel- 
lar medium community, this type of turbulence is called 
super-Alfvenic turbulence, which is favored by some re- 
searchers (see Padoan & Nordlund 1999; Padoan et al. 
2004) as a model of turbulence in molecular clouds. In 
any case, such turbulence is expected to be present in any 
system with magnetic field below the equipartition value, 
which gets subject to intensive driving. As the turbu- 
lence kinetic energy decreases with the scale, we expect the 
Goldreich-Sridhar (1995) model to be valid at some small 
scale when the equipartition is reached, while at larger 
scales we expect hydrodynamic motions to increase the 
energy of magnetic field. 

It is a common knowledge that the effects of mean mag- 
netic field have important astrophysical implications. In 
recent years, interest on the turbulent processes in tan- 
gled magnetic field has been growing. Relevant astrophys- 
ical problems include thermal diffusion in the intraclus- 
ter medium (see Chandran & Cowley 1998; Narayan & 
Medvedev 2001; Cho et al. 2003; Lazarian 2006), cosmic 
ray propagation (see Cassano & Brunetti 2005; Lazarian 
& Beresnyak 2006; Brunetti & Lazarian 2007) as well as 
star formation (Padoan et al. 2004; Li & Nakamura 2004; 
Vazquez-semadeni, Kim & Ballesteros-Paredes 2005). 

Cho & Vishniac (2000a) numerically showed that mag- 
netic energy grows until the magnetic energy density gets 
comparable to the kinetic energy density (see also Kul- 
srud & Anderson 1992; Kulsrud et al. 1997). In this pa- 
per, we present more comprehensive studies on the topic. 
Other aspects of the magnetic field generation also re- 
quire further studies. For example, Haugen & Branden- 
burg (2004) discussed spectral change of MHD turbulence 
by mean field and showed that imposed magnetic field 
lowers the spectral magnetic energy in the inertial range. 
Mac Low (1999) demonstrated that mean magnetic field 
produces anisotropic structures along the mean field di- 
rection in strongly compressible MHD turbulence. Lee 
et al. (2003) discussed the behaviors of energy densities 
and spectral shapes for three different cases (very weak, 
weak, and strong mean magnetic field cases) but only for 
2-dimensional MHD turbulence. They showed the flow 
character in very weak field classes is similar to that of hy- 
drodynamic turbulence, while the strong field cases show 
spectra shallower than the hydrodynamic one. Schekochi- 
hin et al. (2007) discussed the effects of mean fields for 



different values of magnetic Prandtl numbers. In all the 
papers above, the increase of magnetic field energy was no- 
ticed. In this paper, we present a comprehensive study on 
the effect of mean magnetic field on 3-dimensional MHD 
turbulence and turbulence dynamo. 

Another issue that requires clarification is the effect of 
mean magnetic field on the decay of MHD turbulence. 
Stone, Ostriker, & Gammie (1998) and Mac Low et al. 
(1998) numerically showed that damping time-scales of 
compressible MHD turbulence are comparable to the large- 
scale eddy turnover time (see McKee & Ostriker 2007 for 
a collection of related results). Incompressible MHD tur- 
bulence also decays fast (see Cho et al. 2002). While the 
earlier works were mostly focused on two extreme limits - 
zero (see, for example, Biskamp & Muller 2000) and strong 
mean-field limits. In what follows we discuss how the de- 
cay time-scale changes as the mean field strength changes. 

We will first consider the regime of very weak mean field. 
In this regime, we will mainly investigate the growth rate 
of magnetic field. Then, we will consider the effect of in- 
termediate mean fields. We will investigate how magnetic 
and kinetic energy densities, anisotropy, and growth rate of 
magnetic field change with the increase of the strength of 
the mean field. In this paper, we deal with incompressible 
MHD turbulence. We describe our numerical methods in 
§2, present our results for the very weak mean-field regime 
in §3 and discuss the effects of the mean field in §4. We 
compare MHD and hydrodynamic turbulence in §5 and we 
give discussion in §6. We provide our conclusions in §7. 

2. NUMERICAL METHODS 

We have used a pseudospectral code to solve the incom- 
pressible MHD equations in a periodic box of size 27r: 
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= V x (v x B) +?7V 2 B, 
Vv = VB = 0, 



f + VP', (1) 

(2) 
(3) 

where f is random driving force, P' = P + v • v/2, v is 
the velocity, and B is magnetic field divided by {Airp) 1 / 2 . 
We use 21 forcing components with 2 < k < \/l2. Each 
forcing component has correlation time of one. The peak 
of energy injection occurs at k « 2.5. The amplitudes 
of the forcing components are tuned to ensure v ~ 1 for 
the hydrodynamic simulation with v — 0.0074. There- 
fore, one eddy turnover time, ~ L/v, is approximately 2.5 
time units. In this representation, v can be viewed as 
the velocity measured in units of the r.m.s. velocity, v, of 
the system and B as the Alfven speed in the same units. 
Other variables have their usual meaning. The magnetic 
field consists of the uniform background field and the fluc- 
tuating field: B = Bo + b. The Alfven velocity of the 
background field, Bq, varies from to 1. Through out the 

6 The Goldreich & Sridhar model was successfully tested in incompressible 3D MHD simulations in Cho & Vishniac (2000b), Maron & Goldreich 
(2001), Cho, Lazarian & Vishniac (2002), as well as in compressible 3D MHD simulations in Cho & Lazarian (2002, 2003). While the particular 
points of the model, i.e. the particular slope of the spectrum are still the subject of debates (see Maron &; Goldreich 2001; Boldyrev 2005, 2006; 
Beresnyak & Lazarian 2006; Mason, Cattaneo & Boldyrev 2006), the corner stone of the model, which is the critical balance, stays untouched. 

7 In this paper, our goal is to study small-scale turbulence dynamo properties in the asymptotic limit of very small viscosity and very small 
magnetic diffusion. Unit magnetic Prandtl number (y = r]) will be a good approximation in this limit. To support this claim, we ran a 
simulation which is very similar to 256H3 - BolO" 3 but the magnetic Prandtl number (= v/rf) is 0.01. The qualitative behavior of time 
evolution (not shown in this paper) is similar to that of the unit magnetic Prandtl number case: magnetic energy density grows initially and 
saturates later when magnetic energy density becomes comparable to the kinetic one. The average values of v 2 and b 2 are ~ 0.78 and ~ 0.25, 
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paper, we consider only cases where viscosity is equal to 
magnetic diffusion^: 

v = 7]. (4) 

In pseudospectral methods, we calculate the temporal evo- 
lution of the equations |T]) and (J2]) in Fourier space. To 
obtain the Fourier components of nonlinear terms, we first 
calculate them in real space, and transform back into 
Fourier space. We use exactly same forcing terms for all 
simulations. The average kinetic helicity in these simula- 
tions is negative. 

We use an appropriate projection operator to calculate 
VP' term in Fourier space and also to enforce divergence- 
free condition (V • v = V • B = 0). We use up to 384 3 
collocation points. At t = 0, the magnetic field has either 
only uniform component (when Bq ^ 0) or only random 
components (when Bq = 0) and the velocity has a support 
between 2 < k < 4 in wavevector space. 

Either physical viscosity (and diffusion) or hyperviscos- 
ity (and hyperdiffusion) is used for dissipation terms (see 
Table [J). The power of hyperviscosity is set to 3 or 8, such 
that the dissipation term in the above equation is replaced 
with 

-^n(V 2 )"v, (5) 

where n = 3 or 8 and v n is determined from the condition 
v n (N/2) 2n M w 0.5. E Here At is the time step and N 
is the number of grid points in each direction. The same 
expression is used for the magnetic dissipation term. We 
list parameters used for the simulations in Table [TJ We 
use the notation XY-B Z, where X = 384, 320, 256, 144, 
96, or 64 refers to the number of grid points in each spatial 
direction; Y = H or P refers to hyper- or physical viscos- 
ity; Z refers to the strength of the external magnetic fields. 
Diagnostics of our code can be found in Cho & Vishniac 
(2000a). 

We use the following notations: 



B 



- \/ Bq + b 2 : total magnetic field strength or its 
Alfven speed. 



Bq: mean magnetic field or its Alfven speed 
b. 



(= v rms ): the average r.m.s. random magnetic 
field and velocity. Average is taken after turbulence 
reaches a statistically stationary state. 

v^: the zeroth-order magnetic field and velocity 
when Bq = 0. 



E v (k): E v (k) = (l/2)£fc-o.5<fc'<fc+o.5 |v k /| 2 , where 

is the Fourier component of velocity. We define the 
magnetic energy spectrum Eb(k) similarly. 

3. SMALL-SCALE TURBULENCE DYNAMO IN THE VERY 
WEAK MEAN-FIELD LIMIT 

3.1. Expectations 

One of the most important issues in this regime is the 
generation of small-scale random fields from the large-scale 



regular fields. Since large-scale regular magnetic fields are 
observed in almost all astrophysical objects, this issue is of 
great importance in astrophysics. Cho & Vishniac (2000a; 
see also Kulsrud & Anderson 1992; Kulsrud et al. 1997) 
argued that magnetic energy in this regime grows through 
field line stretching and that there are two stages of mag- 
netic field amplification. During the first stage, stretching 
is most active near, or somewhat larger than, the dissipa- 
tion scale (spectral cut-off scale) and the magnetic energy 
spectrum peaks at this scale. It is clear that magnetic en- 
ergy grows exponentially during this stage and that the 
characteristic timescale is the eddy turnover time at the 
dissipation scale. As the magnetic energy grows, the mag- 
netic back reaction becomes important at the dissipation 
scale. When energy equipartition between kinetic energy 
and magnetic energy is reached at this scale, the stretch- 
ing rate slows down and a second stage of slower growth 
begins. Fig. 3 of Cho & Vishniac (2000a) shows that dur- 
ing this stage the peak of the magnetic power spectrum 
moves to larger scales. Using phenomenological arguments 
similar to the ones above, Schekochihin & Cowley (2007) 
argued that magnetic energy grows linearly during the sec- 
ond stage: 

~dT ~ e ' (6) 

where e is the energy injection rate, which should be equal 
to the total energy dissipation rate in a statistically sta- 
tionary state. The linear stage of magnetic energy growth 
ends when stretching on the energy injection scale becomes 
suppressed, which occurs when the magnetic energy den- 
sity becomes comparable to the kinetic energy density. 

3.2. Growth rate at the linear stage 

Let us first consider the growth rate of magnetic en- 
ergy during the linear stage, which has an important con- 
sequence for the strength of the magnetic fields in the 
large-scale structure of the universe (see Ryu et al. 2008). 
Fig. Q] shows time evolution of v 2 and b 2 . All simulations 
started with mean magnetic field strength of Bq = 0.001. 
No random magnetic component was present at the be- 
ginning of the simulation. As simulations go on, random 
magnetic components are generated and magnetic energy 
grows through stretching of magnetic field lines. The 
growth of magnetic energy is slow when viscosity (and 
magnetic diffusivity) is high. For example, Run 64P1- 
_BolO~ 3 shows substantially slower growth rate than Run 
256P-i?olO -3 . The growth rate seems to show a conver- 
gence as viscosity decreases. For example, there is no big 
difference in magnetic field growth rate between Run 256P- 
B Q 10- 3 and Run 256113-BolO -3 . 

We compare the growth rates using simulations with 
different parameters. Right panel of Fig. Q] shows that 
the magnetic energy growth rates during the linear growth 
stage are very similar. Note that we use proper normal- 
ization for both horizontal and vertical axes. We plot only 
high resolution runs. Runs with physical viscosity (and 
magnetic diffusion) show slightly smaller slopes, which is 

BqIO" 3 (see 



respectively. Average is taken over the time interval of (50,110). The magnetic energy density is smaller than that of 256iT3 
Table 1) because magnetic dissipation occurs at smaller wavenumbers. 

8 When a high-order hyperdiffusion is used, the spectral properties near the dissipation cut-off are affected by a strong bottleneck effect. The 
bottleneck effect affects high-fc Fourier components. Since we study mostly the behavior of the total v 2 and b 2 which depend mostly on small-fc 
Fourier components, we believe the bottleneck effect is not a serious issue in our study. 
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reasonable. The strength of mean magnetic field does not 
seem to affect the linear growth rate. However, we can 
clearly observe that, when the mean field is weaker, the 
onset of the linear stage occurs later. It is also worth not- 
ing that even the run with zero mean field (256H8 — BqO) 
shows a similar linear growth rate. In the run with zero 
mean field (256i?8 — BqO), only small-scale (k ~ 70) ran- 
dom magnetic field is present at t=0. 

In code units, the linear growth stage ends at t ~ 40. 
The values of B 2 at that time is ~ 0.4. Therefore, the 
slope during the linear growth stage is around 0.01. But, 
when we represent the slope in terms of normalized energy 
density and time, we obtain different slopes: 

B 2 lt) e t 

~ 0.07— ; —r— + const. 



2E; 



or 



turb 



B 2 (t) 



2E turb (v/L) L/v 



2E, 



0.033 



turb 



L/Vv 2 + b 2 



const., 



(7) 



(8) 



where we use e 
and Eturb = (v 2 



0.16, L ~ 2.5, v ~ 0.9, 



B 2 ~ 1.0, 



B 2 )/2 ~ 0.5 (see right panel of Fig. 1)$ 
A similar linear growth rate has been observed in a recent 
work by Ryu ct al. (2008), in which they derived strength 
of magnetic fields in the large-scale structure of the uni- 
verse. In their model, the linear growth rate derived from 
a simulation plays an essential role. 

3.3. Saturation level 

Using data with relatively low numerical resolutions, 
Cho & Vishniac (2000a) showed that, in the limit of 
v ( = v) ~ > 0, the magnetic energy density in the satu- 
ration stage is comparable to the kinetic energy density. 
This is consistent with the fact that magnetic fields are 
amplified through field line stretching: as we mentioned 
before, stretching of magnetic field at the energy injection 
scale is suppressed only when the magnetic energy den- 
sity becomes comparable to the kinetic energy density. In 
this subsection we present results with higher numerical 
resolutions. 

We list the energy densities in the saturation stage in 
Table 1. We obtained v 2 , b 2 , e, and Dm by averaging over 
(ti, 1%). Here Dm is the magnetic energy dissipation rate. 
It is important to note that these time averages are taken 
after the turbulence has reached a stationary state. Con- 
clusions based on these averaged values do not apply to 
the initial growth phase of the magnetic field. 

Fig. [5] shows v 2 and b 2 as functions of v (= 77) All 
the simulations shown in left panel of Fig. 1 have simi- 
lar kinetic energy densities. However, the magnetic en- 
ergy density obviously depends on the ohmic diffusivity n. 
When v{= rj) goes to zero, the magnetic energy seems to 
approach to ~40% of the total energy. This is somewhat 
larger than the value obtained by Haugen et al. (2003), 
where ~ 30% of the total energy is magnetic. The discrep- 
ancy may stem from the fact that we use hyper-viscosity 
to achieve a very small viscosity (and magnetic diffusion) . 

Runs shown in left panel of Fig. [3] show that the order 
of hyper-diffusion does not strongly affect the growth rate 

9 A careful examination of Fig. 1 reveals that, although v 2 decreases and b 2 increases during the growth stage of magnetic energy, the sum 
of v 2 and b 2 does not change much. In hyper-viscosity runs, the value of v 2 + b 2 is around 1 all the time. This is why we use %/v 2 + b 2 
(a Vv 2 + B 2 ) in eq. (g). 

10 The use of y/u (= y/tf) for horizontal axis is not motivated by theoretical considerations, but by clarity of presentations. 



and saturation level of magnetic energy. Runs shown in 
right panel of Fig.[3]show that numerical resolution slightly 
affects the saturation level of magnetic energy. 

3.4. Exponential growth stage 

Magnetic field is amplified through field line stretching, 
which is initially most active near the dissipation scale. 
As a result, magnetic energy will grow exponentially at 
the beginning. To see this stage more clearly, we plot time 
evolution of magnetic energy in logarithmic scales. We 
use Run 384H3-B 10~ 6 . As expected, Fig. [4] clearly shows 
this exponential growth stage. The strength of the mean 
field is 10~ 6 and the dissipation scale is around k ~ 100. 
When t < 15, the growth rate is exponential and magnetic 
energy spectrum peaks near k ~ 100. At t ~ 15, energy 
equipartition is reached at k ~ 100 and the exponential 
growth stage ends. After t ~ 15, the linear growth stage 
begins and the peak of magnetic energy spectrum moves 
to smaller fc's. 

4. EFFECTS OF THE MEAN FIELD (b ) 

In this section, we consider a similar numerical set-up as 
in the previous section: only a mean magnetic field Bo is 
present at t=0. Then, turbulent motions generate fluctu- 
ating magnetic field b and, as time goes on, the strength 
of the fluctuating field grows. After a certain amount of 
time, which depends on n and B , turbulence reaches a sta- 
tistically stationary state. In this section, we study how 
the strength of the mean field (Bo) affects energy densities 
and other properties of turbulence. As we will see in this 
section, when the mean field is stronger, magnetic energy 
rises faster, the system reaches the statistically stationary 
state more quickly, and magnetic energy density at the 
saturation stage is higher, while kinetic energy density is 
smaller. This result is consistent with an earlier result by 
Tao, Cattaneo, & Vainshtein (1993). Although their main 
conclusion is the suppression of a effect in the presence of 
mean fields, their figures show that magnetic energy den- 
sity at the saturation stage is higher when the mean field 
is stronger. Our result is also consistent with those by 
Schckochihin et al. (2007), where they studied the effects 
of mean fields for different magnetic Prandtl numbers. 



4.1. 



id b 2 



Scaling of energy densities, 1 

Fig. [5] shows time evolution of fluctuating magnetic en- 
ergy density. When the mean field gets stronger, the satu- 
ration level increases and the growth time becomes shorter. 
Runs with both physical and hyper-viscosity show simi- 
lar behavior. In Fig. [5J we explicitly measure the growth 
timescale, which we define the time the energy of the ran- 
dom magnetic field reaches 80% of its saturation value. 
Roughly speaking, Fig. [5] implies that the growth time 
scales with yBo- Numerical resolution or the form of vis- 
cosity is not important unless the mean field is extremely 
small. 

As we can see in Fig. [7J all our simulations, except runs 
with Bo ~ 1, have similar energy dissipation/injection 
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Fig. 1. — Time evolution of kinetic and magnetic energy densities. (Left panel): The level of magnetic energy at the saturation level 
strongly depends on the value of magnetic diffusivity (=viscosity). (Right panel): Comparison of magnetic energy growth rates. The growth 
rates at the linear growth stage are similar. In the case of 256H8-BoO, the mean field is zero and the magnetic energy spectrum at t=0 peaks 
near k ~ 70. In all other cases, only a weak mean field is present at t=0. Note that runs 256H8 — BoO (i.e. run with no imposed field) and 
256H8 — BqW~ 3 show similar growth rates and also similar final saturation levels. 
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Fig. 2. — Normalized average kinetic and magnetic energy densities. Average is taken after turbulence has reached a saturation state. 




Fig. 3. — Comparison of runs. (Left panel): These runs with hyperdiffusion shows that initial strength of the mean magnetic field does not 
strongly affect the growth rate and saturation level of magnetic field. However, amplification of magnetic field (in fact, onset of linear growth 
stage) is delayed when the seed mean magnetic field is weak. The delay will be negligible when the dissipation scale is very small compared 
with the energy injection scale (see i]6.3l for details). (Right panel): These two runs with hyperdiffusion show that numerical resolution also 
does not strongly affect the growth timescale, although two runs show slightly different levels of energy saturation. 



rates e (= f-v). when the mean magnetic field is very weak, the rate is ~ 0.f6 and only weakly depends on v (= r/). 

The rate declines slowly as the external field gets stronger. 
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Fig. 4. — Time evolution of Run 384H3-BolO — 6 . This run with hyperdiffusion shows rapid exponential growth at the beginning, which is 
due to the stretching of magnetic field lines near the dissipation cut-off at kd ~ 100. During this phase, magnetic spectra peaks at k ~ 80 
and moves upwards. After t ~ 15, magnetic energy becomes comparable to the kinetic energy at k^ ~ 100 and the initial rapid exponential 
growth phase ends. 



We list the total energy dissipation rate and the magnetic 
energy dissipation rate in Table 1. In 256P runs, the ra- 
tio -Dm A is ~ 0.7 in weak-mean field cases and it drops 
to ~ 0.55 in strong mean field cases (e.g. 256P--B 0.8 and 
256P-.Boi)- The result for weak mean-field cases is consis- 
tent with Haugen, Brandenburg, & Dobler (2004). 

We measure average energy densities at the statistically 
stationary state and list them in Table 1. In Fig. [51 we plot 
the fluctuating energy densities, v 2 and b 2 (= B 2 — B 2 ), 
as functions of mean field strength Bq. 

First, from the Figure (and Table 1), we note that there 
is no difference between turbulence with no external mag- 
netic field and one with very weak external fields. (For ex- 
ample, compare 256H8-B and 256H8-Bol0" 6 , and 96P- 
B o and 96P-S 10" 3 5 .) When the external fields are 
very weak, the kinetic and magnetic energy densities go 
smoothly to the zero external field limit. From the Fig- 
ure, we also note that the kinetic energy densities are not 
very sensitive to the value of v (= 77) while the magnetic 
energy densities do show a strong dependence on v (= r\) 
(see Cho & Vishniac 2000a for more discussion for Bo sa 
limit). In the case of hyperviscosity simulations, the ratio 
b 2 /v 2 for B -> is greater than 0.6. When B > 0.1, 
magnetic energy, (Bq + b 2 )/2, can be larger than kinetic 
energy, v 2 /2, for 256H3 runs. 

Second, when the external fields are not very strong (say, 
i?o ^ 0.2), the fluctuating energy densities follow 

b 2 cx b^ 2 + c b vB , (9) 

v 2 cx v {0)2 - c v vB , (10) 

where Cb is almost independent of v (= 77) and c v weakly 
depends on v (= ?;) and superscript '(0)' denotes values 
for Bo = 0. Equation ([9]) implies that 

b > {vB ) 1/2 » B . (11) 

4.2. Scaling of total energy 

One consequence of the above scaling relations is that 
the sum v 2 + (c v /cb)b 2 is approximately independent of Bo 
when the external fields are not very strong: 

v 2 + (c v /c b )b 2 ^v^ 2 + (c v /c b )b^ 2 . (12) 



Fig. M shows that total energy (v 2 + b 2 ) does not strongly 
depend on Bo- The constancy of the total energy is es- 
pecially good for runs with high numerical resolutions 
(256P and 256H3 runs). Therefore the Figure implies that 
{cv/cb) ~ 1 for 256P and 256H3 runs. However, the value 
of (c v /cb) is less than 1 for lower resolution runs. Since 
v 2 + b 2 is virtually independent of Bo, we can use the quan- 
tity for normalization. For example, we may define gener- 
alized large-scale eddy turnover time as L/(v 2 + fe 2 ) 1 / 2 . 

4.3. Scaling of energy transfer rate 

Cho and Vishniac (2000a) showed that, when the ex- 
ternal fields are weak/zero, magnetic fields are amplified 
through field line stretching: 

D M oc(v- cB)B 2 /L, (13) 
where Dm is magnetic dissipation and c « 1/1.6. Fig. ITTJ1 
confirms that this result is also true even when the external 
fields are strong. In the stationary state, the magnetic dis- 
sipation (Dm) is balanced by the net energy transferred 
to the magnetic field from the velocity field. The right- 
hand side therefore tells us that the net energy transferred 
to the magnetic field is proportional to large-scale eddy 
turnover rate (v/L) minus an Alfvenic frequency (B / L) 
times a constant. The large-scale eddy turnover rate is 
equal to the stretching rate of the magnetic field when the 
back reaction is zero. We identify the second term on the 
right-hand side of this equation as the effect of the mag- 
netic back reaction. 

4.4. Energy spectra 

We plot energy spectra in Fig.QT] Kinetic spectra peaks 
at k ~ 2.5, which is independent of Bo- Kinetic spectrum 
for the run with Bo = 0.1 (thick solid line) has less energy 
at small k's than that with Bo = 0.001 (thick dashed line). 
On the other hand, magnetic spectrum for B = 0.1 (solid 
line) has more energy at small k's than that for Bo = 0.001 
(dashed line). As the mean field strength grows, the peak 
of the magnetic spectra moves from ~ fc^/2.5 (Bo ~ 
case) to ~ Ul (Bq ) , where k^ is the wavenum- 

ber of the peak of kinetic energy spectra. Kinetic spectra 
are steeper than Kolmogorov spectrum for small k's, while 
magnetic spectra are flatter. Therefore, one should be 
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Fig. 5. — Time evolution of magnetic energy density. The stronger the mean magnetic field, the higher the magnetic energy at the 
saturation stage. The growth time is shorter when the mean magnetic field is stronger. (Left panel): Runs with physical diffusion. (Right 
panel): Runs with hyper-diffusion. 
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Fig. 7. — Average values of the total energy dissipation rate. 



careful when using Kolmogorov spectrum for MHD tur- 
bulence: the kinetic and magnetic spectra are not Kol- 
mogorov for small k's when Bq is weak. Since the slope 
of the total energy spectrum roughly follows Kolmogorov 
slope (see Fig. [TBI and $5]), the slope of kinetic spectra at 
small k's is a function of Bo', when Bq is small, kinetic 



spectra should be steeper than A; -5 / 3 and those for larger 
Bo are closer to fc~ 5 / 3 . When Bo ~ v rm s, the kinetic en- 
ergy spectrum is very close to the Kolmogorov one consis- 
tent with the predictions in Goldreich & Sridhar (1995) as 
well as with some of the earlier simulations (Cho & Vish- 
niac 2000b; Cho et al. 2002). However, we cannot make a 
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stronger conclusion with the present resolution (cf. Maron 
& Goldreich 2001; Beresnyak & Lazarian 2006; Mason et 
al. 2006). 



4.5. Anisotropy 

When B — > v, MHD turbulence tends to be anisotropic 
(Shebalin et al 1983; Goldreich & Sridhar 1995). Here we 



11 In this subsection, wc consider global anisotropy, not local anisotropy. That is, we calculate anisotropy with respect to the direction of the 
imposed field. When we calculate anisotropy with respect to directions that follow the local mean field in smaller sub-volumes, we do obtain 
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Fig. 11. — Energy spectra. Kinetic energy spectrum peaks at the driving scale (=energy injection scale). When the mean field gets stronger 
kinetic spectrum near the driving scale goes down. The magnetic energy spectrum peaks at wavelengths larger than the energy injection scale. 
The location of the peak moves toward smaller wavenumber as the mean field gets stronger. Note that, when Bo < 1, the kinetic spectrum 
does not show a well defined power law and the average slope should be steeper than the Kolmogorov one near the energy injection scale. 



focus on global anisotropy of MHD turbulence with respect 
to the mean magnetic field BoHH Let fc|| be the wavenum- 
ber parallel to the mean field and k± the wavenumber 
perpendicular to the mean field. Matthaeus et al. (1998) 
showed that the anisotropy of MHD turbulence scales 
linearly with the ratio of perturbed and total magnetic 
strength b/B. Their result is mainly for Bo > v. No sys- 
tematic study is available for B < v cases. 

We study anisotropy of MHD turbulence in the B < v 
limit by comparing average k\\ and k±: 



/„*""» d*k fc,||v k | 2 

d*k fc X |v k |2 



(14) 



where k max is the maximum wavenumber and v k is the 
Fourier component of velocity at k. The result in Fig- 
ure [T3 shows that turbulence is virtually isotropic when 
Bo < 0.2v. We note that there is no apparent anisotropy 
of order Bo/ v. When Bo becomes stronger, turbulence 
shows departure from isotropy. Note again that we focus 
on the cases of Bo < v in this paper. It is also noteworthy 
that our result is valid when driving is isotropic. 

4.6. Kinetic helicity 

In Fig. 1131 we plot normalized kinetic helicity, Hk/v 2 , 
as a function of Bo- The kinetic helicity, Hk = (v ■ V x v), 
does show suppression as B increases. But, constancy 
of Hk/v 2 implies that the reduction in Hk is due to the 
reduction in v 2 . Therefore, the figure clearly shows that 
kinetic helicity is not suppressed strongly. Since we can 
write Hk = (v • V x v) ~ k peak ^ Hk v 2 , where k pea k,H h is 
the wavenumber at which kinetic helicity spectrum peaks, 
the constancy of Hk/v 2 means the location of the peak 
wavenumber is not a function of Bo- Indeed the kinetic 
helicity spectra peak at the energy injection scale. 

4.7. Decay timescale 



There have been a lot of studies on the decay law of 
MHD turbulence. Decay law and decay timescale of MHD 
turbulence are of great importance for dynamics of molec- 
ular clouds and star formation. 

Biskamp & Muller (2000) discussed decay of incompress- 
ible MHD turbulence in the limit of zero external field 
cases. However, there has been no systematic studies on 
the decay timescale for various values of the mean field 
strengths. Decay timescale of turbulence energy can be 
estimated by 

(v 2 + b 2 )/e. (15) 

Fig. [33] shows that the rate is not very sensitive to the 
strength of the mean field. Note that, even the case the 
Alfven speed of the mean field (Bo) is comparable to the 
the rms velocity (v ~ 1), the decay rate is not much 
reduced. This fact is consistent with earlier simulations 
of the decaying compressible MHD turbulence (Stone et 
al. 1998; Mac Low et al. 1998) as well as those of incom- 
pressible one (Maron & Goldreich 2001; Cho et al. 2002). 
We present a dimcnsionless dissipation coefficient 

D = (e/2E turb )(L'/u'), ^ (16) 

where e is the energy dissipation rate, E tur b = {v 2 + b 2 )/2, 
3m' 2 = v 2 , and 



L' 



2u' 



k- 1 E v (k)dk. 



(17) 



For 256P MHD runs the value of L' is around ~ 0.76 
and weakly depends on Bo- For example, L' = 0.76 
(256P-B 10^ 3 and 256P-S010" 1 - 5 ), 0.72 (256P-B 10^ 1 ), 
0.68 (256P-B 10-°- 5 ), and 0.71 (256P-B 1). The scale L' 
is a mathematical representation of the energy injection 
scale and the quantity L'/u' is a kind of large scale eddy 
turnover time. Since V ~ 0.76, u' ~ v/V3, for 256P MHD 
runs, we have L'/u 1 ~ 1.5, which is smaller than L/v. In 
256P hydrodynamic run, V ~ 0.62. Note that the defi- 
nition of the dimcnsionless dissipation rate is not exactly 
the same as that in Kaneda et al. (2003), in which they 
used 



{e/u"){L'/u'). 



(18) 



anisotropy for all values of Bq. This kind of anisotropy is local anisotropy. Scale-dependent anisotropy discussed in Goldreich & Sridhar (1995) 
is local anisotropy (see Cho & Vishniac 2000b and Cho et al. 2002 for details). Local anisotropy is important for local physical processes 
(e.g. pitch-angle scattering). But, what matters for external observers is mostly global anisotropy (or global isotropy). 
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The dimensionless dissipation coefficient is known to 
be around 0.5 — 0.6 both for hydrodynamic turbulence 
(Kaneda et al. 2003) and MHD turbulence (Stone et 
al. 1998; Haugen, Brandenburg, & Dobler 2004; see McKee 
& Ostriker 2007). Our results are consistent with earlier 
findings: the value is around 0.5 — 0.6 when the mean field 
is weak. This value somewhat drops when mean field is 
stronger. This has to do with the fact that the Kolmogorov 
constant depends on the value of the mean field (see next 
section). 

5. COMPARISON WITH HYDRODYNAMIC SIMULATIONS 

It is often claimed that MHD turbulence with weak im- 
posed magnetic field is similar to hydrodynamic turbu- 
lence. Is this true? When we introduce a mean field in a 
fully turbulent medium, what will happen to energy densi- 
ties? In this section, we compare hydrodynamic and MHD 
turbulence. 

We first compare energy densities. We run hydrody- 
namic and MHD simulations with the same initial condi- 
tions, except Bq, and the same prescribed driving force. 
Then we compare energy densities at the saturation stage. 
Fig. fTS] shows the behavior of the energy densities as a 
function of v (= rf). It is interesting that < v > in hy- 
drodynamic turbulence and < v 2 + 1.6b 2 > in MHD cases 
scale similarly. Right panel of the figure shows that the 
similarity is also true for individual time basis. 

The claim has a long history that total energy spectrum 
in MHD turbulence follows a Kolmogorov spectrum (see, 
for example, Kida, Yanase, & Mizushima 1991; Miiller & 
Biskamp 2000). We plot the energy spectra of hydrody- 
namic and MHD turbulence in Fig. [16] The compensated 
spectra in right panel of the figure shows that the total en- 
ergy does show a slope compatible with the Kolmogorov 
one. However, the compensated spectrum of MHD turbu- 
lence can be higher or lower than that of hydrodynamic 
one depending on the mean field strength. In fact, the 
vertical location in the right panel of Fig. [TBI corresponds 
to the Kolmogorov constant Ck ■ 

E(k) = C K e 2/3 k- 5/3 . (19) 
The Kolmogorov constant for hydrodynamic turbulence is 
around 1.6 (see, for example, Yeung & Zhou 1996). On 
the other hand, in MHD case, Biskamp & Miiller (2000) 
reported that Ck ~ 2.3 in their decaying turbulence sim- 
ulations. Our results show that the constant can be a 
function of Bq. 

We compare the dimensionless dissipation coefficient 
(equation |16| ) for hydrodynamic and MHD turbulence in 
Fig. [171 We define the coefficient as in equation (fT6j) . For 
MHD cases we use E tur b = {v 2 + b 2 )/2 and, for hydrody- 
namic cases, Eturb = v 2 /2. In the plot, we also mark the 
coefficient defined as in equation (fl8|) . where u' 2 is used 
instead of 2E tU rb, for hydrodynamic cases (see squares on 
the Figure). In hydrodynamic runs, the coefficient is de- 
pendent on the value of v. However, in MHD runs, the 
coefficient does not show strong dependence on v (— rj). 

6. DISCUSSION 

6.1. Magnetic reconnection 

Astrophysical magnetic fields are really ubiquitous. The 
origin of the large-scale regular magnetic fields is uncer- 
tain. Nevertheless, the turbulence dynamo discussed in 



this paper provides a simple and efficient way to produce 
small-scale random fields. It only requires the existence of 
turbulence, which is as ubiquitous as magnetic fields them- 
selves (see Armstrong, Spangler & Rickett 1995; Lazarian 
& Pogosyan 2000; McKee & Ostriker 2007). Nevertheless, 
small-scale turbulence dynamo assumes fast magnetic re- 
connection. 

The reconnection is, indeed, fast in the numerical sim- 
ulation that we deal with here. If, in astrophysical fluids, 
magnetic reconnection is slow, i.e. happens on time scales 
longer than the dynamical time, both the models of strong 
turbulence and the generation of magnetic field by turbu- 
lence are misrepresented by the numerics. In this situation 
we expect the formation of unresolved magnetic knots that 
should substantially alter the dynamics of the fluid and 
magnetic fields. If, however, reconnection if fast, i.e. hap- 
pens on the time scales comparable to the Alfvenic wave 
propagation time scale, then our calculations do represent 
generation of astrophysical magnetic fields. We have am- 
ple astrophysical evidence and some models of fast recon- 
nection are rather robust. For instance, a model of turbu- 
lent reconnection in Lazarian & Vishniac (1999; see also 
Vishniac, Lazarian & Cho 2003; Lazarian, Vishniac & Cho 
2004) predicts that magnetic field change its topology over 
the eddy turnover time. If, however, magnetic reconnec- 
tion is not a robust process, i.e. is fast only in particular, 
e.g. collisionless environments (see Shay et al. 1999), the 
turbulent generation of magnetic field will differ substan- 
tially from one environment to another. Future testing of 
magnetic reconnection should clarify this issue. 

6.2. Turbulence with moderate imposed fields in 
astrophysical flows 

We have considered the effects of weak mean magnetic 
fields (0 < B /^/i¥p < v rms ) on MHD turbulence. Note 
that in our simulations < v 2 >— v 2 ms ~ 1 and effectively 
we can assume p = 1 for incompressible runs. In our code 
units, the factor 47r disappears. Therefore, the parameter 
space we have considered is < Bq < 1. This work pro- 
vides missing link between two extreme limits: Bq — > 
limit and Bq / \JAnp > v rrns limit. 

The limit of Bq — > has been studied by many re- 
searchers (Pouquet & Patterson 1978; Meneguzzi, Frisch, 
& Pouquet 1981; Kida, Yanase, & Mizushima 1991; Cho 
& Vishniac 2000a). In this paper, we have shown that 
there is no difference between turbulence with no external 
magnetic field and one with very weak external fields. In 
this sense, the limit of Bq — > is not a special limit. Cho 
& Vishniac (2000a) showed that energy equipartition be- 
tween small scale kinetic energy and total magnetic energy 
occurs at a scale about 3 times smaller than the energy 
injection scale. Magnetic spectrum of hyperviscosity sim- 
ulation peaks at a wavelength about 2-3 time larger the 
kinetic energy peak. These results imply that, when we 
have a tangled magnetic field with no mean magnetic field 
in a fully developed turbulence, the characteristic size, Is, 
of the magnetic field is about 2-3 times smaller than the 
energy injection scale of the turbulence. 

When < Bo/yJAirp < v rms , the turbulence is called 
super- Alfvenic in the interstellar medium community. The 
intracluster medium may fall in this regime (Ryu et 
al. 2008). Although it is still uncertain, the interstellar 
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Fig. 14. — (Left panel): Decay timescale. (Left panel): Dimcnsionless dissipation rate. Note that the inverse of the dimcnsionlcss dissipation 
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Fig. 15. — Comparison between MHD and hydrodynamic turbulence. (Left panel): Average energy density at the saturation stage. In 
MHD cases, Bo = 0.001 for all simulations. The runs that correspond to v K, arc 256H3-Hydro and 256H3-BolO~ 3 . (Right panel): The 
similarity between v? dro and (v 2 + 1.66 2 ) mhd holds true even for time evolution. 




Fig. 16. — Comparison of energy spectra. (Left panel): Total energy spectra for MHD and kinetic energy spectrum for hydro- 
dynamic turbulence. Right panel shows that the slope is compatible with Kolmogorov spectrum. (Right panel): Compensated spec- 
tra show that, when Bo is weak, the total energy spectrum of MHD is lower than hydrodynamic one. This is understandable because 
(v 2 + o 2 )mhd < (v 2 + 1.&} 2 )mhd ~ v hydro regardless of the strength of the mean magnetic field and because e's in MHD with Bo ~ 
and hydrodynamic cases are nearly same. However, when Bo ~ v, e reduces and, therefore, the compensated spectrum is higher than the 
hydrodynamic one. We plot spectra averaged over time. 



medium and/or molecular clouds in Galaxy may also fall in gave an estimate of field strength derived from many dif- 
this regime (Padoan et al. 2004a; Beck 2001). Beck (2001) ferent observations. He proposed that the regular com- 
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Fig. 17. — Dimensionless energy dissipation coefficient. Both hydrodynamic and MHD cases are shown. In MHD cases, So = 0.001 for all 
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ponent in our galaxy is ~ 4^G and the total component 
~ 6/iG. On the other hand, Heiles & Crutcher (2007) sug- 
gested that the energy ratio Et U rb/E m ag for cold neutral 
medium (CNM) is ~ 1.3 in our galaxy. Fig. [5] shows 
that, when Bq > 0.1, B 2 can be larger than v 2 for 256H3 
runs. Therefore, we expect that the mean field in the in- 
terstellar medium should be larger than ~ 0.1. Note that, 
when Bo is larger than ~ 0.2, turbulence becomes glob- 
ally anisotropic (Fig. [T2"]) . Therefore, ISM turbulence is 
either marginally anisotropic if Bq ~ 0.1 or anisotropic if 
Bq > 0.2. Turbulence in this case will reach saturation 
level very quickly within 3 or 4 eddy turnover times at 
most. The correlation length scale of magnetic field will 
be very close to the energy injection scale. 

6.3. Model for growth of magnetic field in astrophysical 

flows 

In numerical simulations, the dissipation scale cannot be 
arbitrarily small. Instead, the dissipation scale is limited 
by the numerical resolution. In actual astrophysical fluids, 
the dissipation scale will be much smaller than the one we 
can achieve with numerical simulations. Then what will 
be the time evolution of the magnetic energy in an actual 
astrophysical fluid which has a very high Reynolds num- 
ber? Here, we only consider the case that Bo is extremely 
small. When numerical resolution or Reynolds number is 
high, the dissipation scale is small and the eddy turnover 
time at the dissipation scale is short. As a result, the ex- 
ponential growth is fast. Suppose that the linear growth 
stage begins at t = t'. We plot the situation in left panel 
of Fig. [TBI Time evolution of magnetic energy for t > t' is 

B 2 (t)-v 2 d = C(t-t'), (20) 

where, C is the slope at the linear growth stage. The ki- 
netic energy at the dissipation scale, ~ v d , and the time t' 
are given by 

v 2 ~v 2 (l d /L) 2 /*~ v 2 (R e )-V 2 , (21) 

and 

t' ~T d ln(v d /B ) 2 
~ teddy(ld/L) 2/3 HRe- 1/4 v/B ) 2 
~ teddy (Re)- 1 ' 2 HRe-^v/Bo) 2 , (22) 



where t ec id y is the large scale eddy turnover time and 
we used the fact IdVd/v ~ 1 ~ (Lv/v)(ld/L)(vd/v) ~ 
(Re) (I4/L) 4 ' 3 . Here Id is the dissipation scale, v the 
large scale rms velocity, v d the velocity at the dissipation 
scale, Td is proportional to the eddy turnover time at the 
dissipation scale, and Re = Lvjv the Reynolds number. 
Eqs. IpTj) and ([22]) imply that both v 2 d and t' go to zero 
when Reynolds number is very large. Therefore, we get 

B 2 (t) » Ct (23) 
in real astrophysical fluids with very large Reynolds num- 
bers. When numerical resolution is finite, time evolution 
of B 2 will look like 

B 2 (t)KC(t-t') (24) 

because v 2 d is negligible for reasonably high numerical res- 
olutions. 

Results in £)3.4I (see left panel of Fig. [3J) show that the 
onset of the linear stage occurs later when the mean field 
is weak. Then how large is the delay? We plot the situa- 
tion in right panel of Fig. [18] Note that the linear stage 
begins when energy equipartition is reached at the dissipa- 
tion scale, which should be independent of the strength of 
the mean field. When the mean field is strong, the linear 
stage begins at t\. When the mean field is weak, the linear 
stage begins later at i 2 . Roughly speaking, the magnetic 
energy at the exponential growth stage is 

B 2 (t)^B 2 cxp(-t/T d ). (25) 
Therefore, we have 

B 2 )X cxp(-i 1 /T d ) = B 2 2 exp(-t 2 /r d ). (26) 
Solving this equation, we get 

t 2 - ti cx Tdln(B 0) i/So )2 ) 

cx t eddy (Re)- 1 ' 2 ln(S a/Bo, 2 ). (27) 
Of course, when the Reynolds number is large, i 2 — t\ is 
nearly zero, which mean that the strength of the mean 
field does not matter much in astrophysical fluids. 

7. CONCLUSION 

We have studied the growth of magnetic field in the 
presence of weak/moderate mean magnetic fields. We 
have considered only the cases of v — rj. We have found 
that equipartition between the kinetic and magnetic en- 
ergy densities occurs at a scale somewhat smaller than the 
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Fig. 18. — Growth of magnetic field energy in the presence of very weak mean magnetic field. (Left panel): Initially magnetic energy grows 
exponentially. E-folding time in this stage is the eddy turnover time at the dissipation scale. When the magnetic energy reaches equipartition 
with kinetic energy at the dissipation scale, exponential stage ends and linear growth stage begins. In actual astrophysical fluids with high 
Reynolds number, t' and are very small and, therefore, we expect B 2 (t) oc t for most of the time until the saturation stage. (Right panel): 
When the mean field strengths are different, the linear growth stage begins at different times. When the Reynolds number is high, ti — t\ 
will be very small. 



kinetic energy peak, which is consistent with previous re- 
sults with lower numerical resolution by Cho & Vishniac 
(2000a). We have found that runs with different numeri- 
cal resolutions and simulation parameters show consistent 
results for the slope of the linear growth stage. In real as- 
trophysical flows with very weak mean field and very large 
Reynolds number, the growth of magnetic energy will fol- 
low 

V ' - 0.033—^===. (28) 



2E, 



turb 



L/Vv 2 + B 2 



We have also studied the effects of external magnetic 
fields on MHD turbulence. We have shown that mag- 
netic energy density increases and kinetic energy density 
decreases as the external field becomes stronger. To be 
specific, we have shown that, when Bq < 0.2, 

b 2 - b^ 2 cx vB Q , (29) 

v 2 - w (0)2 cx -vB , (30) 

where b 2 is the magnetic energy density in the presence of 
the external field and b^ 2 is the value when the external 
field is zero. 

When the external magnetic field is not very strong 
(~ less than 0.2 times the rms velocity), the turbulence 



remains statistically isotropic, i.e. there is no apparent 
anisotropy of order Bq/v. The slope of kinetic spectra at 
small k's is a function of Bq: when Bq is small, kinetic 
spectra should be steeper than fc -5 / 3 and those for larger 
Bq are similar to fc -5 / 3 . 

We have compared MHD and hydrodynamic turbulence. 
Spectra of total energy in MHD turbulence are similar to 
those of hydrodynamic turbulence. The Kolmogorov con- 
stant in MHD cases can be higher or lower than that of 
hydrodynamic turbulence, depending on the strength of 
the mean field. 
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